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When the interbrane separation in compact Randall-Sundrum models is stabilized via the 
Goldberger-Wise mechanism, a potential is generated for the four dimensional field that encodes this 
geometric information, the so-called radion. Due to its origin as a part of the full five dimensional 
metric, the radion couples directly to particles on both branes. We exhibit the exponential growth 
in the number of brane particles due to parametric amplification from radion oscillations and discuss 
some of the consequences of this process for brane cosmology. 
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Compact Randall-Sundrum models (RSI) provide 
an extremely interesting relation between hierarchies in 
four dimensional physics and the warping of the five di- 
mensional bulk geometry in which our universe is em- 
bedded. These models consist of two 3-branes placed 
at the fixed points of an S ,1 /Z 2 orbifold, which serve as 
the boundaries of a slice of a five dimensional anti-de 
Sitter space-time (AdSs). In order to balance the bulk 
cosmological constant, the branes must have equal and 
opposite tensions, which are in turn related to the value 
of the bulk cosmological constant. 

These models necessarily contain a geometrical mode 
that encodes the physical separation between the branes, 
the so-called radion. We can write the full 5-d line ele- 
ment in a way that exhibits it explicitly O , 

ds 2 = e -Hy+f(*y 2ky ]g^( X )dx»dx» 

+ {l + 2kf(x)e 2k y) 2 dy 2 . (1) 

The radion is given by the function f(x) and the metric 
is written in such a way as to properly incorporate the 
symmetries of the geometric setup. In particular, this 
form guarantees that the radion will not mix with mass- 
less degrees of freedom such as the 3 + 1 gravitons Q . 
From the linearized Einstein equations 

Rnu - \g^R = o v 2 / = o, (2) 

where R^ v and are the Ricci tensor and covariant 
derivative associated with the metric g^(x), we can see 
that the radion behaves as a massless 4-d field. This is 
in accord with its status as the Goldstone boson corre- 
sponding to the fact that the action is invariant under 
the operation of placing the second brane at any dis- 
tance from the first. Furthermore, it should be noted 
that the radion has no Kaluza-Klein tower associated 
with it @,|| . 

The Goldstone nature of the radion precludes the sys- 
tem from picking a value for the interbrane separation. 
Physically, this separation controls the hierarchy between 
particle physics scales and that of gravity, so that in or- 
der to make the RSI scenario work, some means of fixing 



the expectation value of the radion must be given. This 
is the essence of the Goldberger-Wise (GW) mechanism 
H which introduces a bulk scalar with a non-trivial pro- 
file in the bulk and potentials on each of the branes. For 
simplicity, the self-couplings of the GW scalar were cho- 
sen large enough so as to "stiffen" the potentials to the 
point that the value of the GW scalar on each brane 
is completely fixed. Integrating the action of the GW 
scalar over the bulk coordinate gives rise to a potential 
for the radion, whose expectation value can then be fixed 
in terms of the parameters of the GW action to take the 
correct value so as to enforce the required hierarchy. 

Our interest in the radion was piqued by the fact that 
since it is a mode contained in the full metric, it should 
couple to all brane fields (or at least those whose actions 
break Weyl invariance on the brane) in a universal man- 
ner. From a cosmological perspective, we should expect 
that the radion will not settle to its minimum instanta- 
neously, but that it will oscillate for some time about the 
minimum of the potential generated by the GW mech- 
anism. As with the inflaton at the end of the inflation- 
ary era H, these oscillations can induce the explosive 
production of brane particles via parametric amplifica- 
tion, i.e. preheating Q. The energy density and equation 
of state of these particles will influence the cosmological 
evolution of the relevant brane. We calculate this particle 
production below for the case of a brane scalar coupled 
to the radion. 

We expect that when the back reaction on the radion 
due to this particle production is taken into account, the 
radion evolution will become dissipative and its expec- 
tation value will settle to its minimum. However, this 
is beyond the scope of this paper and will be treated, 
together with other related topics, in Q. 

In order to calculate brane preheating effects due to 
radion oscillations, we first need the radion action. The 
kinetic terms come from the five dimensional Einstein 
action, supplemented by the Gibbons-Hawking terms on 
the boundary branes which contribute to the effective po- 
tential for the radion. The GW action is that of a free 
scalar field of mass m with potentials to fix the value of 
the field to be vqM^ 2 , and v\M^ 2 on the UV and the 
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IR branes, respectively, where M5 is the five dimensional 
Planck mass. The calculation of the radion effective ac- 
tion obtained by integrating out the extra dimension is 
relatively straightforward and will be given in detail in 
; the resulting effective Lagrangian is 
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We have defined the following quantities in Eq. (0): 
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where Ay is the location of the IR brane relative to the 
UV brane at y = 0, A w is the cutoff for the effective 
theory of the radion such that at energies above Aw the 
radion dynamics become strong and r\ (< 1) parameter- 
izes the change in the tension of the IR brane in order 
to account for the effects of the GW scalar. Rcf. ||] also 
found an O(TeV) mass for the radion based on a calcu- 
lation in the full theory rather than the effective action 
approach which yielded Eq. (||). 

We envisage a scenario in which at the early stages of 
an RS universe, when the typical temperature of the bulk 
is still above e~ kAy M 5 , the space-time is generally unsta- 
ble to the formation of a black hole horizon in the bulk 
whose Hawking temperature is the same as that of the 
bulk. The model then resembles the Randall-Sundrum 
model [|| with a single UV brane and a horizon. As the 
universe cools below the energy scale associated with the 
IR brane, the free energy of the model with two branes 
becomes lower JhJ. After its formation it will generally 
not be at its equilibrium position but will oscillate about 
it with a natural amplitude of the order of Ayy, thus 
driving the brane preheating process. 

Thus consider a scalar field Q(x^) confined to the IR 
brane at y = Ay. Its action involves the radion mode 
through the induced metric. In the small amplitude 
regime for $, we can write the scalar action as 



5$ = / d x ■ 
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where h a b is the induced metric on the IR brane. Using 
Eq. (pi), we can rewrite this as 
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where we have performed the rescalings $ = e~ kAy <& 
and m$ = e~ kAy rh<s> so as to obtain a canonically nor- 
malized kinetic term when the radion is at its equilibrium 
value F = 0. The couplings of the radion to Standard 
Model fields were also considered in Q where precision 
electroweak observables were studied. In contrast, we 
are more concerned with dynamical aspects of these cou- 
plings of the radion to brane fields. 

The oscillating radion acts as a time dependent back- 
ground in which $ evolves, i.e. we neglect radion quan- 
tum fluctuations as well as the effects of the back-reaction 
due to the produced particles or other fields. This should 
be a reasonable assumption when the back-reaction on 
the radion from other fields is small and the higher or- 
der terms in the effective radion Lagrangian Eq. (^) can 
be neglected. The time dependence of the radion back- 
ground and its couplings to <I> imply that energy can be 
drained off the radion to produce 4> quanta. In the ab- 
sence of back-reaction this process proceeds indefinitely. 
A more realistic expectation is that the back-reaction will 
damp out the motion of the radion (7). We take the ra- 
dion background to be of the form 



F(t) = F A W cos(mF^), 



(7) 



where F is dimensionless. 

To compute the number of particles produced, we ex- 
pand $ in terms of creation and annihilation operators, 
at and for a mode with momentum fc, 
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The mode functions arc then determined by the equation 
of motion for the scalar field from Eq. (Q), 

d 2 <&- d<I>- 

-^-± + F o m F sm(m F t)^+u;l(t)<f> ; : = 0, (9) 
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where Lo\{t) = \k\ 2 + 

We shall track the number N^(t) of quanta of the initial 
state defined via 



%(0) = -«^(0)%(0). 



(10) 



This defines a vacuum state |0) which is annihilated by 
the operators a^. 

The periodicity of the coefficient functions in Eq. (||) 
under t — > t + 2nmp 1 implies, by Floquet theory (see 
e.g. |nj) that there will be values of the dimensionless 
momentum |fc|/m$ for which <I>^ (i) will undergo expo- 
nential growth, i.e. there will be unstable bands. We 
display this in Fig. |l|. 

For momenta in the unstable bands, the particle num- 
ber grows exponentially, as can be seen from Fig. |^, where 
we plot the expression Q for N^(t), 

m) = 
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FIG. 1: The broad instability band for the differential equa- 
tion Eq. (|^) when Fa = 0.5. The dark region shows where the 
amplitude of &e{t) grows exponentially. Both tuf and k = \k\ 
are expressed in units of m$. 
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FIG. 2: lniVg(t) for m F = 4m* and F = 0.5. Note that for 
momenta within the instability band of Fig. the number of 
particles produced grows exponentially; t and k are expressed 
in units of mj and m$, respectively. 
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This explosive particle production cannot help but af- 
fect brane physics. We shall explore the cosmological 
effects of such particle production in detail in H , but we 
can make some observations at this point. First, these 
are standard model particles presumably being produced 
starting when the brane temperature is O(TeV). While 
they begin with a non-thermal equation of state, we 
should expect that by the time the universe cools down 
to nucleosynthesis temperatures, they will have thermal- 
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FIG. 3: The equation of state w = p/p as a function of time. 
The dotted line indicates average value of w at later times. 
For this case, m_p = 4m$, this average is approximately 0.25. 



t is given in units of m^ 1 . 



ized. However, depending on the abundance of these par- 
ticles, as well as their equation of state, they can modify 
the expansion law of the IR brane, i.e. our universe. To 
follow this effect carefully would require computing the 
expectation value of the stress tensor of the produced 
particles and using it to source the relevant FRW-likc 
equations that describe the expansion of the brane. We 
have estimated the density and pressure of the produced 
particles and plotted the behavior of the "equation of 
state" w = p/p as a function of time in Fig. ||. 

The initial radion amplitude will fix the location of 
the instability band as will the ratio tof/to$. This will 
change the (time averaged) equation of state accordingly 
and thus the expansion rate. 

The initial distribution of the Standard Model fields on 
the brane depends on the details of the mechanism that 
produces the IR brane. If the Standard Model fields are 
essentially in their vacuum state, then the subsequent 
evolution of the non-thermal particles could differ sub- 
stantially from the standard picture of a radiation dom- 
inated universe with a temperature < O(TeV). Even 
with an initial thermal population of fields, the paramet- 
ric amplification due to the radion could still produce a 
significant fraction of the total energy density in a non- 
equilibrium distribution. Such a large non-equilibrium 
component at electroweak energies could produce inter- 
esting effects. 

We have seen that radion oscillations, corresponding 
to a "breather mode" for the two branes, can give rise to 
explosive particle production on the IR brane. Presum- 
ably, similar effects could arise on the UV brane as well, 
although the effective field theory techniques used here 
may not apply. Preheating on the UV brane may have 
even greater consequences for cosmology on the IR brane, 
since energy densities on the UV brane can have an ex- 
ponentially large effect on IR brane cosmology |p| . Ra- 
dion induced preheating could also be important in elec- 
troweak physics, and in particular in electroweak baryo- 
genesis [H, where this particle production could affect 
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sphaeleron physics. We leave all this to our larger work 
§• 
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